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ABSTRACT 

We present results from the first three-dimensional calculations ever to follow the 
collapse of a molecular cloud core (~ 1CT 18 g cm -3 ) to stellar densities (> 0.01 g cm -3 ). 
The calculations resolve structures over 7 orders of magnitude in spatial extent (~ 
5000 AU — 0.1 R©), and over 17 orders of magnitude in density contrast. With these 
calculations, we consider whether fragmentation to form a close binary stellar system 
can occur during the second collapse phase. We find that, if the quasistatic core that 
forms before the second collapse phase is dynamically unstable to the growth of non- 
axisymmetric perturbations, the angular momentum extracted from the central regions 
of the core, via gravitational torques, is sufficient to prevent fragmentation and the 
formation of a close binary during the subsequent second collapse. 

Subject headings: accretion - hydrodynamics - ISM: clouds - methods: numerical - 
stars: formation 
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1. Introduction 

The formation of close binary stellar systems is an, 
as yet, unsolved problem in the field of star forma- 
tion. Calculations following the isothermal collapse 
and fragmentation of molecular cloud cores can suc- 
cessfully produce binaries with separations <;l AU 
(e.g. Boss 1986). However, once a collapsing molecu- 
lar cloud core becomes optically thick, the gas tem- 
perature increases with density, a hydrostatic core is 
formed with radius ~ 4 AU ( Larson 196S ), and frag- 
mentation is suppressed (Boss 1986, 1988). 

A potential opportunity to form close binaries via 
fragmentation occurs when the temperature in the 
hydrostatic core reaches ~ 2000 K and molecular 
hydrogen begins to dissociate, resulting in a second 
collapse within the hydrostatic core (Larson 1969). 
When the dissociation is complete, the collapse is fi- 
nally stopped with the formation of a second hydro- 
static, or stellar, core. Fragmentation during the sec- 
ond collapse has been obtained by Boss (1989) and 
Bonnell & Bate (1994), although in the former paper 
the binary quickly merged. In both papers, fragmen- 
tation was found to occur only if the ratio of the rota- 
tional energy to the magnitude the gravitational po- 
tential energy, /?, exceeded the value required for dy- 
namic instability to non-axisymmetric perturbations 
(J3 ^ 0.274; Durisen et al. 1986) after a stellar core 
formed. Otherwise, only a stellar core and stable disc 
form. However, in both papers, only the inner regions 
of the first hydrostatic core were modelled and the ini- 
tial conditions were chosen somewhat arbitrarily. Ide- 
ally, to study the formation of close binary systems, 
one would wish to follow the collapse of an optically- 
thin, molecular cloud core all the way to stellar den- 
sities. However, such a calculation requires resolving 
densities from 10~ 18 — 0.01 g cm -3 and length scales 
from ~ 10 17 — 10 10 cm in three dimensions. 

To achieve such high spatial resolution in three 
dimensions we use the Lagrangian smoothed parti- 
cle hydrodynamics (SPH) method. With SPH, the 
spatial resolution is given by the smoothing length 
(over which hydrodynamical properties are calcu- 
lated) which in modern implementations is continu- 
ously variable (e.g. Benz et al. 1990). The smoothing 
lengths are set by ensuring that each particle contains 
a certain number of neighbours, N ne igb-> or equiva- 
lently a fixed mass, within two smoothing lengths. 
Thus, in contrast to a grid-based code which has 
spatially-limited resolution, SPH has mass-limited 
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Fig. 1. — The gas temperature T (solid line), Jeans 
mass Mj (dotted line), and Jeans radius Rj (dashed 
line) , as a function of density in a collapsing molecular 
cloud core (adapted from Tohline 1982). The mini- 
mum Jeans mass of « 4 x 10~ 4 M Q occurs at the end 
of the second collapse phase when the dissociation of 
molecular hydrogen is almost completed at a density 
of ~ 10~ 3 g cm~ 3 . 

resolution which automatically gives greater spatial 
resolution in regions of higher density. Bate & Burk- 
ert (1997) recently demonstrated this mass-limited 
resolution. They found that to follow the collapse and 
fragmentation of a Jeans-unstable gas cloud reliably, 
the minimum resolvable mass should always be less 
than the Jeans mass (see also Whitworth 1998). In 
practice, this means that a Jeans mass should always 
be represented by at least ~ 27V ne i g h particles. 

The mass-limited resolution of SPH is ideal for 
studying the collapse and fragmentation of molecular 
cloud cores because there is a minimum Jeans mass 
in the problem (Figure 1). By contrast, there is no 
minimum Jeans length. This is a problem for grid- 
based codes which must resort to nested or adaptive 
grids (e.g. Burkert & Bodenheimer 1993; Truelove et 
al. 1998). With SPH, if the number of particles used is 
sufficient to resolve the minimum Jeans mass, a calcu- 
lation can be followed to arbitrary densities with the 
required spatial resolution given automatically with 
increasing density. 

2. Method 

The calculations were performed using a three- 
dimensional SPH c ode based on a version orig inally 
developed by Benz (|Benz 1990| , |Benz et al. 1990|) . The 
code uses a tree to calculate the gravitational forces 
and find nearest neighbours. The smoothing lengths 
are variable in time and space, subject to the con- 
straint that the number of neighbours of each particle 
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must remain approximately constant at iV nc i g h = 50. 
The gravitational force between particles is softened 
over the hydrodynamic smoothing length using spline 
softening derived from the smoothing kernel. Hence, 
the gravitational softening length and hydrodynamic 
smoothing length are equal and decrease with increas- 
ing density. The standard form of the artificial viscos- 
ity is used ( Monaghan fc Gingold 1983| ), with the pa- 
rameters a v = 0.5 and f3 v = 1.0. The SPH equations 
are integrated using a second-order Runge-Kutta- 
Fehlberg integrat or. Individual time-steps ar e used 
for each particle ( Bate, Bonncll fc Price 1995 ). Fur- 
ther details can be found in the above references. 

The code does not include radiative transfer. In- 
stead, to model the behaviour of the gas during the 
different phases of collapse, we use a piece- wise poly- 
tropic equation of state P = Kp 7 , where P is the 
pressure, p is the density, K gives the entropy of the 
gas, and the ratio of specific heats, 7, is varied as 
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where the densities are in g cm -3 (see Figure 1). The 
values of 7 and the transition densities are derived 
from Tohline (1982). The variable value of 7 mim- 
ics the following behaviour of the gas. The collapse 
is isothermal (7=1) until the gas becomes optically 
thick to infrared radiation at p w 10~ 13 g cm~ 3 , be- 
yond which 7 = 7/5 (appropriate for a diatomic gas). 
When the gas reaches a temperature of « 2000 K 
(p = 5.7 x 10~ 8 g cm -3 ), molecular hydrogen be- 
gins to dissociate and the temperature only slowly 
increases with density. In this phase we use 7 = 1.15 
to model both the decreasing mean molecular weight 
and the slow increase of temperature with density, 
the latter of which has an effective 7 ps 1.10. Finally, 
when the gas has fully dissociated (p ~ 10~ 3 g cm -3 ), 
the gas is monatomic and 7 = 5/3. The value of K is 
defined such that when the gas is isothermal, K — 
with c s = 2.0 x 10 4 cm s -1 , and when 7 changes the 
pressure is continuous. 

To test that the above equation of state captures 
the important elements of the gas's behaviour, we 
have performed spherically-symmetric, one-dimensional 
(1-D), finite-difference calculations of the collapse of 
a molecular cloud core using the above equation of 
state and obtained excellent agreement with results 
from 1-D calculations incorporating radiative trans- 
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Fig. 2. — Maximum density (solid line) and maximum 
temperature (dotted line) versus time for the collaps- 
ing molecular cloud core. Time is given in units of the 
initial free-fall time (tg = 1.785 x 10 12 s). The right 
graph has expanded axes to show the second collapse 
phase in greater detail. 

fer (e.g. Larson 1969; Winkler & Newman 1980a, b). 
We have also tested that the three-dimensional (3-D) 
SPH code can indeed accurately resolve the collapse 
down to stellar densities by performing comparison 
calculations with the same initial conditions as used 
by the 1-D finite-difference calculations. So long as 
the number of particles used by an SPH calculation 
is sufficient to resolve the Jeans mass throughout the 
calculation (see below) , we obtain excellent agreement 
with the results from the 1-D finite-difference code. 
The details of these test calculations will be presented 
in a later paper. 

3. Calculations 

We present results from a 3-D SPH calculation 
following the collapse of an initially uniform-density, 
uniform-rotating, molecular cloud core of mass M = 
1 Mq and radius R = 7 x 10 16 cm. The angular fre- 
quency is f2 — 7.6 x 10 -14 rad s _1 , and the ratios of 
the thermal and rotational energies to the magnitude 
of the gravitational potential energy are a = 0.54 
and (3 = 0.005, respectively. To satisfy the resolu- 
tion criterion of Bate & Burkert (1997), we require 
that the minimum Jeans mass during the calculation 
(« 4 X 10~ 4 M Q ) contains at least « 2iV ncigh = 100 
particles. Hence, we use 3 x 10 5 equal-mass particles. 

The evolution of the calculation is as follows (Fig- 
ure 2). The initial collapse is isothermal. When the 
density surpasses 10" 13 g cm" 3 , the gas in the cen- 
ter is assumed to become optically thick to infrared 
radiation and begins to heat (t = 1.009 tg). The 
heating stops the collapse at the center and the first 
hydrostatic core is formed (t = 1.015 tg) with maxi- 
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mum density « 2 x 10~ n g cm -3 , mass w O.Ol M Q 
(« 3 x 10 3 particles), and radius « 7 AU. As the 
first core accretes, its maximum density only slowly 
increases at first. However, the first core is rapidly 
rotating, oblate, and has « 0.34 > 0.274, mak- 
ing it dynamically unstable to the growth of non- 
axisymmetric perturbations (e.g. Durisen et al. 1986). 
At t « 1.023 tff, after about 3 rotations, the first 
core becomes violently bar-unstable and forms trail- 
ing spiral arms. This leads to a rapid increase in 
maximum density as angular momentum is removed 
from the central regions of the first core (now a 
disc with spiral structure) by gravitational torques 
(t > 1.015 tff). When the maximum temperature 
reaches 2000 K, molecular hydrogen begins to disso- 
ciate, resulting in a rapid second collapse to stellar 
densities (t = 1.030 tff). The collapse is again halted 
at a density of w 0.007 g cm -3 with the formation 
of the second hydrostatic, or stellar, core. The initial 
mass and radius of the stellar core are « 0.0015 M 
(« 5 x 10 2 particles) and « 0.8 R©, respectively. 
Finally, an inner circumstellar disc begins to form 
around the stellar object, within the region under- 
going second collapse. The calculation is stopped 
when the stellar object has a mass of « 0.004 M© 
(w 1.2 x 10 3 particles), the inner circumstellar disc 
has extended out to « 0.1 AU, and the outer disc 
(the remnant of the first hydrostatic core) contains 
w 0.08 Mq (« 2.4 x 10 4 particles) and extends out 
to « 60 AU. Note that the massive outer disc forms 
before the stellar core. This final state is illustrated 
in Figures 3 and 4. If the calculation was evolved 
further, the inner circumstellar disc would continue 
to grow in radius as gas with higher angular momen- 
tum fell in. Eventually, the inner disc would meet the 
outer disc with only a small molecular dissociation 
region between the two. 

The value of during the second collapse, or after 
the stellar core and disc have formed, never exceeds 
the value necessary for dynamic instability to non- 
axisymmetric perturbations (0 = 0.274). Therefore, 
unlike in the calculations of Bonnell & Bate (1994), 
fragmentation to form a close binary system cannot 
occur. A similar result is obtained with SI = 3.4 x 
10~ 14 rad s- 1 and = 0.001 initially. 

Since whether or not fragmentation can occur de- 
pends critically on the value of 0, it is important that 
there has not been significant spurious transport of 
angular momentum by artificial viscosity or the diffu- 
sion of SPH particles. As SPH is a Lagrangian tech- 



nique, it is relatively simple to keep track of the grav- 
itational, pressure, and viscous torques on each parti- 
cle. Thus, we can determine the cause of the change 
in angular momentum of a volume of gas between any 
two times. As one might expect, we find that there is 
virtually no angular momentum transport during dy- 
namic collapse. Viscous evolution is more likely dur- 
ing the quasistatic evolution of the first core between 
its formation and the second collapse. In Figure 5, we 
consider the cause of the angular momentum trans- 
port during the quasistatic first core phase. Viscous 
evolution, which has presumably been active during 
the « 3 rotations that the first core made before it 
became violently bar-unstable, is hardly detectable. 
By contrast, it takes only 1/4 of a rotation for grav- 
itational torques to remove most of the angular mo- 
mentum from the gas in the center of the core. 

The absence of fragmentation could also be due to 
the value of 7 during the second collapse being over- 
estimated. This would lead to the value of beta that 
is achieved at the end of the collapse being underes- 
timated. To determine the dependence of the results 
on 7, we performed calculations with 7 as low as 1.05. 
This corresponds to an isothermal collapse but with 
the mean molecular weight still decreasing due to the 
dissociation of molecular hydrogen. In this extreme 
case does briefly exceed 0.274 as the stellar core be- 
gins to form, but drops below the critical value again 
over approximately one rotation as the core increases 
in mass and before non-axisymmetric perturbations 
are able to grow significantly. The inner disc that 
forms around the stellar core also has a brief period 
with > 0.274 until more mass is accreted. Weak, 
transient spiral arms are visible in the disc over sev- 
eral rotations, however, there is no tendency for frag- 
mentation. 

4. Conclusions 

Larson (1969) performed the first hydrodynamical 
calculations of the collapse of a molecular cloud core 
to stellar densities using a spherically-symmetric, one- 
dimensional code. After nearly 30 years, it is now pos- 
sible to perform three-dimensional calculations. Us- 
ing SPH, we have followed the three-dimensional col- 
lapse of a molecular cloud core over 17 orders of mag- 
nitude in density contrast and 7 orders of magnitude 
in spatial extent, resolving both the isothermal and 
second collapse phases, and the formation of the first 
and second hydrostatic cores. 
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Using this technique, we have begun a study of 
whether or not close binary stellar systems may be 
formed during the second collapse phase of a molecu- 
lar cloud core, as proposed by Bonnell & Bate (1994). 
We find that, if the first hydrostatic core is formed 
with sufficient rotational energy to make it dynami- 
cally unstable to the growth of non-axisymmetric per- 
turbations, fragmentation during the second collapse 
is prevented because of the removal of angular mo- 
mentum from the central regions of the first core via 
gravitational torques. The formation of close binaries 
during the second collapse phase may be possible for 
molecular cloud cores with different initial conditions, 
especially if they are rotating slowly enough that the 
first core is rotationally stable. Different initial con- 
ditions will be investigated in a subsequent paper. 
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Fig. 3. — The state of the system at the end of the calculation. The six panels on the left give the density and 
velocity in the plane perpendicular to the rotation axis and through the stellar core. The six panels on the right 
give the density and velocity in a section down the rotation axis. In each case, the six consecutive panels give the 
structure on a spatial scale that is 10 times smaller than the previous panel to resolve structure from 3000 AU to 
i=a 0.2 R Q . The remnant of the first hydrostatic core (now a disc with spiral structure), the inner circumstellar disc, 
and the stellar core are all clearly visible. The distance scale (in AU), velocity scale (in km s _1 ), and logarithm 
of the minimum density (in g cm -3 ) are given under each panel. The maximum density is 0.03 g cm~ 3 and the 
logarithm of density is plotted with contours every 0.5 dex. 
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Fig. 4. — The state of the system at the end of the calculation. The left graph gives the density (solid line) and 
smoothing length (dashed line) as a function of radius and the mass enclosed within radius r of the center of the 
stellar core (dotted line). The right graph gives the radial (solid line) and rotational (dotted line) velocities as 
functions of radius from the center of the stellar core. The densities, smoothing lengths and velocities are the mean 
values in the plane perpendicular to the rotation axis and through the stellar core. The stellar core (r < 0.004 
AU), inner circumstellar disc (0.004 < r < 0.1 AU), region undergoing the second collapse (0.1 < r < 1 AU), outer 
circumstellar disc formed from the first core (1 < r < 60 AU), and isothermal collapse region (r > 60 AU) are 
clearly visible. 
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Fig. 5. — The cause of the change in angular momentum during the calculation. We show the system at t = 
1.0234 tg. The first core has undergone approximately 3 rotations since its formation; it has been violently bar- 
unstable for just 1/4 of a rotation. The left graph gives the density as a function of radius (solid line) and the 
mass enclosed within radius r of the center of the first core (dotted line). The right graph shows the fractional 
change, since the beginning of the calculation, of the angular momentum of the gas as a function of radius (solid 
line), and the contributions to this change from gravitational (short-dashed line), pressure (dotted line) and viscous 
(long-dashed line) torques. Gravitational torques from the bar instability dominate the loss of angular momentum 
from the center of the first core; this angular momentum is deposited in the outer parts of the core. Pressure forces, 
due to the spiral arms, act to oppose the angular momentum transport. Viscous transport is negligible. 



7 



